Fundamental Uncertainty in the BAO Scale from Isocurvature Modes 
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Small fractions of isocurvature perturbations correlated with the dominant adiabatic mode are 
shown to be a significant primordial systematic for future Baryon Acoustic Oscillation (BAO) sur- 
veys, distorting the standard ruler distance by broadening and shifting the peak in the galaxy 
correlation function. Untreated this systematic leads to biases that can exceed 10a in the dark 
energy parameters even for Planck- level isocurvature constraints. Accounting for the isocurvature 
modes corrects for this bias but degrades the dark energy figure of merit by at least 50%. The 
BAO data in turn provides extremely powerful new constraints on the nature of the primordial 
perturbations. Future large galaxy surveys will thus be powerful probes of the earliest phase of the 
universe in addition to helping pin down the nature of dark energy. 



Introduction Large galaxy surveys fall squarely in the 
realm of astrophysics and traditionally we think of them 
as living almost independently of the physics of the early 
universe. In this Letter we show that this assumption 
breaks down rather dramatically in the case of Baryon 
Acoustic Oscillations (BAO) (for a recent review see |l| 
). BAO surveys are a key component of the global plan 
for the next two decades in cosmology because they are 
believed to provide a robust and powerful statistical stan- 
dard ruler that can probe dark energy. They have shown 
to be robust to a variety of potential systematic effects 
which only become important at the ~ 1% level @, 
However, here we show that there is a much more signifi- 
cant "systematic" arising from the possibility of isocurva- 
ture modes correlated with the dominant adiabatic per- 
turbation which may have been generated during the 
early universe. 

To understand this systematic, consider the standard 
ruler provided by the distance that sound waves can prop- 
agate in the primordial plasma. The standard picture 
based on adiabatic perturbations suggests that the key 
scale is the sound horizon: 



c s (l + z)dt 



Cs(z') 

H(z'Y 



dz , 



where c s (z) = l/y/3 (1 + R b {l + z)" 1 ) and R b = 
31500o;b (T cm h/2.7K) -4 . The measurement of the late- 
time clustering of galaxies in the transverse direction 
probes the angular diameter distance given by d^{z) — 
rj_(l + z)~ 1 /AQ where r± is the intrinsic size of r s in 
the transverse direction and AO is the position of the 
peak in the angular correlation function, while the clus- 
tering on a scale ru along the line of sight probes the 
Hubble parameter, H(z) = Az/rn. Measurements of the 
angular diameter distance and Hubble parameter in a 



series of redshift bins using the BAO technique provide 
an effective probe of the properties of dark energy 0, 0] , 
with prospective constraints on the equation of state of 
dark energy, wo, and its evolution w a , as low as 0.01 and 
0.05, respectively, for a future space-based spectroscopic 
mission (e.g., Adept [5|), with forecasts for current ex- 
periments (e.g., Boss at the level of 0.03 and 0.07 
respectively. 

Systematic effects that affect the position and shape 
of the Baryon Acoustic Peak (BAP), such as nonlinearity 
and redshift-space distortions, have been studied and can 
be treated without a significant impact on dark energy 
constraints 0, 0] • Here we concentrate on the possibil- 
ity that the initial conditions were not purely adiabatic. 
Isocurvature modes deform this characteristic scale which 
manifests in the anisotropies in the cosmic microwave 
background and the clustering of matter. We investigate 
whether a small admixture of isocurvature modes could 
potentially degrade constraints on dark energy parame- 
ters if allowed for, or bias the measured values if not taken 
into account. Put more generally, we investigate whether 
it is possible to decouple the physics of the generation of 
the primordial density perturbation from the constraints 
arising from our observations at late times, even with 
strong prior constraints on the isocurvature modes from 
CMB data. Constraints from Wmap 3 year data indicate 
that a 50% admixture of three isocurvature modes with 
the adiabatic mode is permitted Q , whereas forecasts for 
the Planck experiment indicate that isocurvature mode 
admixtures will be constrained to below the 10% level [§]. 

What are the possible origins of isocurvature modes? 
The simplest possibility perhaps is multiple field inflation 
3 . Ilfj| . with the curvaton mechanism as a special case 
111 ]. The resulting isocurvature perturbation is a lead- 
ing candidate to explain any primordial non-Gaussianity 
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and can, in certain cases, explain the observed asymme- 
try in the CMB [12]. While the simplest, adiabatic mod- 
els of inflation are currently preferred [13| , it is possible 
that some isocurvature contamination will be uncovered 
in future experiments and indeed this would be very for- 
tuitous since it would provide new handles on the physics 
of the very early universe. In this letter we argue that al- 
lowing for the possibility of isocurvature modes is crucial 
in future BAO surveys and that as a reward, such surveys 
can provide a powerful lens on the early universe. 
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FIG. 1: Baryon mass profile for adiabatic and isocurvature 
modes at decoupling, with the correct relative amplitude be- 
tween the AD (black solid), NID (blue short dot-dashed line) 
NIV (green long dot-dashed line), BI (cyan dashed line) and 
CI (red dotted line) modes. 
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FIG. 2: Baryon mass profile for adiabatic and isocurvature 
modes at z — 0, with the same labelling as in figure 1. For 
clarity the NID, NIV and CI mass profiles have been multi- 
plied by a factor of two, while the BI mass profile has been 
multiplied by a factor of six. 

The BAO peak with adiabatic and isocurvature initial 
conditions The features of the BAO peak are not only 
sensitive to the background dynamics, but also to the 



evolution of the cosmic perturbations, in particular the 
manner in which the initial relative perturbations be- 
tween the different species were established. In addition 
to adiabatic fluctuations, in which the equation of state 
of the universe is spatially constant and the curvature 
is perturbed, there is the possibility of regular isocur- 
vature perturbations, in which the equation of state be- 
tween the different species varies to keep the curvature 
constant, and of correlations between the adiabatic and 
isocurvature perturbations 



Acoustic oscillations in the photon-baryon fluid are de- 
scribed by the photon density evolution equation in har- 
monic space 



R X , 2 2 r ^ 

— S 7 + k c s 5 7 = -- 



R 



1 + R 



h + h 



= F(k,r), 



where R = 



Rb 



is the baryon-to-photon density ra- 



( 1 + z ) 

tio, h is the metric field in synchronous gauge, and the 
dot refers to the conformal time derivative. Solving this 
equation prior to decoupling we find that 

S 7 (k, r) — As sin kr s (r) + Ac cos kr s {f) 
+A T f (l + i?(r')) 1/2 sin[fcr s (r)-fcr s (T')]F(fc,T')dr', 



with the initial conditions that define the regular modes 
given by 
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where Sl C; o and f^o are, respectively, the cold dark mat- 
ter and the baryon densities today, with R v and R 1 the 
fractional energy densities of neutrinos and photons at 
early times, respectively. The NIV mode starts with a 
non-zero perturbation in the photon velocity so that it 
stimulates the sm(kr s ) harmonic in <5 7 , in contrast to the 
NID mode which has an initial non-zero perturbation in 
the photon density and excites the cos(kr s ) harmonic. 
We note that the adiabatic mode is sourced purely by 
the gravitational driving term F(k,r), which is constant 
on large scales and at early times, so that the adiabatic 
solution excites a cos(fcr s ) harmonic, similar to the NID 
mode. At later times, though, this approximation breaks 
down when the mode enters the horizon, and the adi- 
abatic solution deviates from a pure cosine mode, with 
the consequence that its acoustic peaks are offset in phase 
from the NID peaks. In the transition to matter dom- 
ination, the gravitational driving term switches on as a 
source for the CI and BI modes, whereas its contribution 
to the NID and NIV modes remains negligible. 
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The different harmonics stimulated by the various 
modes result in a relative shift of the BAO peak posi- 
tion. We observe the same shift in the series of acous- 
tic peaks exhibited in the cosmic microwave background 
anisotropy spectrum , though in this case it is recom- 
bination, as opposed to decoupling, that sets the acous- 
tic scale. To quantify the shifts in the BAO peak of the 
isocurvature modes relative to the adiabatic mode, in fig- 
ure Q] we plot the baryon mass profile []J§] of the various 
modes at decoupling, defined by 



M b (r,z)= / T b (k,z) 



sin kr 
kr 



k\2 e -k 2 a 2 /2 dk 



where Tf, = 6b/ k 2 is the baryon transfer function. Prior 
to decoupling and on large scales, 6b is exactly |<5 7 for all 
modes except the BI and NID modes in which the baryon 
and photon density perturbations differ by a constant, 
respectively, 1 and i?„/i? 7 . 

The mass profile captures the evolution of an initial 
point-like perturbation at the origin, though in practice 
we use a narrow gaussian, with width tr _1 in k space. 
The initial density perturbation expands out as a spher- 
ical wave at the sound speed [ijj so that at decoupling 
there is an excess of baryons at the sound horizon scale, 
r s w c s Tdcc- Note that even though the photons decou- 
ple from the baryons at z « 1080, it is only at z ~ 500 
that the baryons stall, due to the fact that the growing 
mode is dominated by the velocity field on small scales, 
which does not decay instantaneously after decoupling. 
This sets the scale for the low-redshift BAO peak, which 
is shown in figure El We observe that the BAO feature 
for the isocurvature modes differs in shape and position 
from the adiabatic BAO peak. In the case of the NID 
and NIV modes there is a pronounced peak which is off- 
set from the adiabatic BAO peak due to the coupling to 
different harmonics, as described above. In the case of 
the BI and CI modes the acoustic wave has merely im- 
printed a ripple onto the homogeneous sea of baryons at 
decoupling, which evolves into a knee in the baryon mass 
profile at the BAO scale at late times. It is also interest- 
ing to note that the amplitude of the mass perturbation 
at the origin is much smaller for the isocurvature modes 
because the curvature, and thus the mass fluctuation, is 
initially unperturbed. It is clear that small admixtures of 
the isocurvature modes and their cross-correlations can 
distort the shape and location of the adiabatic BAO peak. 

Dark energy constraints Assuming that the variations 
in the cosmological parameters (including those model- 
ing isocurvature) are small, we can model the likelihood 
function of a dataset as a multivariate Gaussian centered 
on a fiducial adiabatic ACDM Universe. Based on the 
noise estimates for the BOSS and Planck experiments, 
we can compute estimates of the errors on the cosmolog- 
ical parameters using a Fisher matrix formalism, by per- 
turbing the cosmology around the fiducial model. When 
perturbing the dark energy model away from A, we allow 



0.7 



-0.1 



-0.9 











ADEPT bios 














AD: 


Plonck+BOSS 












_ ISO: 


Plonck+BOSS 












'. AD:' 


Plonck+ADEPT 












_ ISO: 


Plonck+ADEPT 




























— 1 .00 —0.97 


-0.94 








0.20 
















0.10 




























BOSS 


bias 


















0.00 
















-0.10 
















-0.20 











1.20 



1 .03 



-0.87 



-0.70 



FIG. 3: la error ellipses for (wo,w a ) from the Fisher ma- 
trix calculation for the Boss (blue) and Adept (red) experi- 
ments with the Planck data as a prior, assuming adiabaticity 
(solid) and for the fully correlated isocurvature case (dashed). 
The fiducial model is marked by a black star. The average 
biases in the dark energy parameters that could potentially 
incurred as a result of the incorrect assumption of adiaticity 
are shown for the two experiments by the dotted lines. All 
other parameters have been marginalized over. An enlarged 
version of the error ellipses is shown in the insert. 



for dynamics and parameterize its equation of state us- 
ing w(a) — wq + (1 — a)w a where a = 1/(1 + z) HI, 19j - 



To model deviations away from adiabaticity, we adopt 
the isocurvature parameterization implemented in [20j, 
where the different modes and their cross-correlations are 
described by 9 parameters, measuring the fractional con- 
tributions of the various correlations (auto and cross) to 
the overall total power spectrum. Note that for the com- 
putation of the Planck Fisher matrix, we follow [2l| and 
re-introduce the strict geometric degeneracy between the 
dark energy density Qx and wq, w a which may be artifi- 
cially broken in the standard Fisher matrix computation, 
leading to under-estimates of errors. 

We are concerned with how well the equation of state 
parameters wo and w a can be measured. In order to 
quantify the constraining power of the data, we com- 
pute the Dark Energy Task Force (DETF) figure of merit 
(FoM), which is defined as the reciprocal of the area in 
the wq — w a plane, enclosing the 95% confidence limit 
(CL) region 2l|. We are concerned with the change in 
the FoM when isocurvature is introduced relative to the 
case of pure adiabaticity. 

We compute the potential errors on wq and w a for the 
case of pure adiabaticity and for the scenario where all 
isocurvature modes and their cross-correlations are ad- 
mitted, while marginalizing over all other cosmological 
parameters. The results for both the Boss and Adept 
experiments are shown in the first and last rows of table 
HI The Boss FoM is found to increase by 50% when this 
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Experiment Boss Adept 

Modes wo w a wo w a 

AD 0.023 0.071 0.0068 0.028 

AD+CI+ ( AD, CI) 0.024 0.071 0.0069 0.030 

AD+NID+(AD,7VLD) 0.024 0.072 0.0069 0.028 

KD+mV+{AD,NIV) 0.025 0.072 0.0069 0.031 

ISO (ALL) 0.031 0.087 0.0075 0.037 



TABLE I: Table summarizing the constraints on (wo, w a ) for 
adiabatic and an admixtures of uncorrelated adiabatic and 
isocurvature modes, marginalizing other all other parameters, 
for the Boss and Adept experiments. The fiducial model 
assumes adiabaticity. 

additional freedom is introduced in the initial conditions, 
while the Adept FoM degrades by 90%. These changes 
are illustrated in figure [3] which compares the 68% con- 
fidence regions for the adiabatic and the fully correlated 
isocurvature cases. Evidently, no single mode and its cor- 
relation are responsible for the change in the allowable 
(u>o, w a ) region, but rather a mixture of all extra degrees 
of freedom. Note that the forecasted errors on wq and 
w a for the adiabatic case in our analysis are considerably 
smaller than those reported in likely due to our as- 
sumption of spatial flatness, but what is important is the 
relative change in the FoM. 

The question that we wish to ask is what bias in the 
estimates of the dark energy parameters could potentially 
be induced by incorrectly assuming adiabaticity? For a 
Gaussian-distributed likelihood function, it can be shown 
that the linear bias in a set of parameters that we wish 
to constrain, 50i, due to erroneous values of a set of fixed 
parameters, 6(f>j, is [22| 

where F is the Fisher sub-matrix for the parameters 
we wish to constrain and F 9 ^ is a Fisher sub-matrix 
constructed from the product of the derivatives of the 
power spectrum with respect to the parameters being 
constrained and those which are being fixed. In our 
case j labels the nine isocurvature mode amplitudes, in- 
correctly fixed to zero, k labels the eight cosmological 
parameters that are biased, and i labels the subset of 
two dark energy parameters whose bias is of interest to 
us. In order to set 5(f>j, we diagonalize the combined 
Planck and large-scale structure (LSS) Fisher matrix 
and select the eigenvector, ei with the smallest eigen- 
value Xi . This corresponds to the direction in parameter 
space which is least constrained by the data. We then 

take 8(f> j = y^f^e,;. For the BOSS experiment, we find 
that Swo = -0.044 and Sw a = 0.12, while for Adept the 
biases are found to be Swq = —0.0095 and 5w a — —0.061. 

However this is only one particular direction that 
weakly constrains all the parameters, not necessarily the 



dark energy parameters. In order to explore the full 
range of the bias, we use a set of 10,000 random lin- 
ear combinations of the eigenvectors to compute 5(f) j and 
the corresponding biases. We find that the mean bi- 
ases are fi(Swo) = —0.13 and fi(8w a ) — —0.87 with 
<t(5wo) = 0.072 and a(6w a ) = 0.39. The implication is 
that if the initial conditions of our Universe are comprised 
of a sub-dominant contribution from isocurvature modes 
(within the la constraints from Planck and the selected 
LSS survey) , the assumption of adiabaticity could lead to 
an incorrect 6a detection of non-A dark energy model or 
a 12cr false claim of dynamics. Alternatively, A could 
be found to be consistent with the data when in fact 
w(z) 7^—1. The potential bias incurred by the adiabatic 
assumption in the case of a more advanced BAO exper- 
iment such as Adept is somewhat smaller for wq, with 
H(5w ) = 0.021 and a(Sw ) = 0.0072, while the mea- 
surement of w a could be inaccurate at the level of 17a 
with fi(Sw a ) = 0.48 and a(8w a ) — 0.18. The reduction in 
the bias is encouraging as one would expect that with an 
increase in the constraining power of the survey comes 
a higher risk of making false claims. We note that the 
estimates in this paper so far are optimistic because we 
assume perfect knowledge of the galaxy bias and redshift 
distortions. 

Although allowing for isocurvature modes degrades the 
dark energy constraints relative to the perfectly adia- 
batic case, there is a powerful positive. The volume of 
the 9-dimensional isocurvature Fisher ellipse is roughly 
2— 4x 10 9 smaller than that from Planck alone, showing 
that using CMB plus BAO data in union provides excep- 
tionally good constraints on the early universe relative to 
the CMB alone. More specifically, we find that the error 
bars on the isocurvature parameters decrease by 30% to 
as much as 100% for certain modes when the LSS data 
(either Boss or Adept) is added to the Planck data. 

Conclusions With forecasted constraints on dark en- 
ergy from BAO experiments at the level of a few percent 
made possible by the large volumes probed by the most 
recent generation of redshift surveys, it is important to 
explore the full spectrum of possible BAO systematics. 
In this Letter, we have revisited the assumption of adia- 
batic initial conditions. We have found that the admis- 
sion of isocurvature admixtures alters the position of the 
BAP, and by assuming adiabaticity we run the risk of 
incorrectly attributing a shift in the peak (away from the 
predicted value in a ACDM model) to the presence of 
dark energy. 

Ignoring isocurvature modes can substantially bias the 
estimates of the dark energy parameters, leading to a 6a 
(3a) incorrect measurement of wo or as much as a 12a 
(17a) bias in w a for BOSS (Adept) . Allowing for general 
initial conditions removes this bias at the expense of an 
increase in the 2a region in (w , w a ) space by 50% (90%) 
for the BOSS (Adept) survey indicating that the as- 
sumption of adiabaticity can lead to an under-estimation 
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of the errors on the dark energy parameters. On the other 
hand, BAO data substantially improves the constraints 
on isocurvature modes. 

Acknowledgements CZ is funded by a NRF/DST 
(SA) Innovation Fellowship and a National Science Foun- 
dation (USA) fellowship under grant PIRE-0507768. PO 
is funded by the SKA (SA). SMK received support via the 
Meraka Institute via funding for the South African Cen- 
tre for High Performance Computing (CHPC). BB and 
KM acknowledge support from the National Research 
Foundation, SA. 



[1] B.A. Bassett & R. Hlozek, larXiv:0910. 52241 (2010) 
[2] H-J. Seo & D. Eisenstein, Ap. J., 598, 720 (2003) 
[3] H-J. Seo & D. Eisenstein, Ap. J., 665, 14 (2007) 
[4] C. Blake & K. Glazebrook, Ap.J. 594, 665 (2003) 

L. Amendola et al, Mon. Not. Roy. Astron. Soc. 357 

(2005) 

A. Rassat et al, Phys. Rev. D 80, 123516 (2009) 
[5] E. Sefusatti & E. Komatsu, Phys. Rev. D76, 083004 



[6 
[7; 

[8 

[9 
[10 

[11 
[12 

[13 

[14 
[15 
[16 

[17; 

[is: 

[19 
[20 
[21 
[22 



(2007) 

D. Schlegel, et al. larXiv:0902.4680l (2009) 

R. Bean et al, Phys. Rev. D74, 063503 (2006) R. Bean et 

al, Phys. Rev. D 74, 063503 (2006) 

M. Bucher et al, Phys. Rev. Lett. 87, 191301 (2001) 

S. Mollerach et al., Phys. Rev. D 44, 1670 (1991) 

D. Polarski & A. A. Starobinsky, Phys. Rev. D 50, 6123- 

6129 (1994) 

D.H. Lyth & D. Wands, Phys. Lett. B 524, 5-14 (2002) 
A.L. Erickcek et al, Phys. Rev. D 80, 083507 (2009) 
J. Valiviita & T. Giannantonio, Phys. Rev. D 80, 123516 
(2009) 

M. Bucher et al, Phys. Rev. D 62, 083508 (2000) 
J. Dunkley et al., Phys. Rev. Lett. 95, 261303 (2005) 

D. Eisenstein, et al, Ap. J., 664, 660 (2007) 

S. Bashinsky & E. Berschinger, Phys. Rev. D65, 123008 
(2002) 

M. Chevallier & D. Polarski, Int. J. Mod. Phys. D 10, 
213 (2001) 

E. V. Linder, Phys. Rev. Lett, 90, 091301 (2003) 
K. Moodley, et al. Phys. Rev. D 70, 10 3520 (2004) 
A. Albrecht et al, |arXiv:astro-ph/060959T1 (2006) 

A.N Taylor et al. , Mon. Not. Roy. Astron. Soc. 374, 
1377 (2007) 



